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Abstract: 

In this paper we report numerical and experimental results on the scaling prop- 
erties of the velocity turbulent fields in several fiows. The limits of a new form of 
scaling, named Extended Self Similarity (ESS), are discussed. We show that, when a 
mean shear is absent, the self scaling exponents are universal and they do not depend 
on the specific flow (3D homogeneous turbulence, thermal convection , MHD). In 
contrast, ESS is not observed when a strong shear is present. We propose a general- 
ized version of self scaling which extends down to the smallest resolvable scales even 
in cases where ESS is not present. This new scaling is checked in several laboratory 
and numerical experiment. A possible theoretical interpretation is also proposed. 
A synthetic turbulent signal having most of the properties of a real one has been 
generated. 

Submitted for publication in Physica D September 1995 
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1 Introduction 



In order to characterize the statistical properties of fully developed turbulence 
one usually studies the scaling properties of moments of velocity differences at the 
scale r: 

Sp{r) =< \v{x + r) - v{x)f >=< \5v{r)f > (1) 

where < ■ ■ ■ > stands for ensemble average and v is the velocity component parallel 
to r. At high Reynolds number Re = UqL/v the Sp{r) satisfies the relation 

Sp{r) oc r^(P) (2) 

for L > r >> r/fc where L is the integral scale, rj^ = (z/^/e)^/"^ is the dissipative 
(Kolmogorov) scale, e is the mean energy dissipation rate, u the kinematic viscosity 
and Uq the R.M.S. velocity of the flow. The range of length L > r >> rj^, where the 
scaling relation (0) is observed, is called the inertial range. The Kolmogorov (K41) 
theory predicts ({p) = p/3, but experimental]^ and numerical |Q results show 
that ({p) deviates substantially from the linear law. This phenomenon is believed to 
be produced by the intermittent behaviour of the energy dissipation [^] which can 
be taken into account by rewriting eq.(0) in the following way: 

Sp{r) oc< ePj^ > rP/^ oc t^^p/^'^+p/^ (3) 

where is the average of the local energy dissipation e(x) on a volume of size r 
centered on a point x. A comparison of eq.(|l|) and eq.(||) leads to the conclusion 
that the scaling exponents r(p/3) of the energy dissipation are related to those of 
SphYCip)=rip/3)+p/3. 

Since the Kolmogorov (K62) theory Q many other models, ||^, 0, |^, [ID 



TT[| have been suggested to describe the behaviour of the C(p)- However, it turns 
out that the ({p) may be not universal in non homogeneous, anisotropic flow and 
may depend on the location where measurements are done. Specifically, they may 
have different values if one measures either far away from boundaries, where tur- 
bulence is almost homogeneous and isotropic, or in locations of the flow where a 
strong mean shear is present. The ({p) depend also on the way in which turbulence 
is produced, for example 3D homogeneous turbulence, boundary layer turbulence, 
thermal convection and MHD. Thus there is the fundamental question of under- 
standing in which way all these parameters influence the scaling laws. Furthermore 
all the above mentioned models assume the existence of two well defined intervals 
of lengths that are the inertial range and a dissipation range. According to idea 
of multiscaling these two ranges may eventually be connected by an intermediate 
region where the viscosity begins to act |jl2|. However this idea of a well defined 
inertial range, where viscosity does not act at all, and the idea of multiscaling turns 
out to be incompatible with the recently introduced new form of scaling, which has 
been named Extended Self Similarity (ESS) [|T^ [|l4|(see section 2 below). 
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ESS has been observed in 3D homogeneous and isotropic turbulence both at low 
and high Re and for a wide range of scales r with respect to scaling . In contrast 
ESS is not observed when a strong mean shear is present ||15[. All these experimental 
observations show that also the mechanisms by which energy is actually dissipated 
in a flow are very poorly understood. Specifically one would like to understand how 
viscosity acts on different scales. This is clearly an important point in order to safely 
use large eddy simulations in real applications. 

The purpose of this paper is to rationalize all the above mentioned results on 
scaling both in presence and absence of a shear. We propose a generalized form of 
ESS which has been checked in many different flows. We have also generated a signal 
which has all the statistical properties of a real turbulent signal. Our interpretation 
of ESS and of this generalized scaling suggest that there is no sharp viscous cut-off 
in the intermittent transfer of energy. 

The paper is organized as follows: in section 2 we remind the properties of ESS, 
in section 3 we discuss the systems where ESS is not observed, in section 4 the 
hierarchy of structure functions is described, in section 5 the generalized form of 
scaling is discussed, in section 6 a possible theoretical interpretation is proposed, in 
section 7 we discuss the multiscaling. Finally conclusions are given in section 8. 



2 Extended Self Similarity 

Extended Self Similarity (ESS) is a property of velocity structure functions of ho- 



mogeneous and isotropic turbulence |13, 14|. It has been shown using experimental 



and numerical data |]T6| that the structure functions present an extended scaling 



range when one plots one structure function against the other, namely: 

Sn{r) oc 5^(r)^("''-) (4) 



where (3{n,m) = Q{n) / C,{m) . The details of ESS have been reported elsewhere |jT4|. 
In the following we describe only the main features. 

As an example we consider two experimental data sets at different R\, which 
is the Reynolds number based on the Taylor scale {R\ ~ 1.4 Re^^"^) The two 
experiments are a jet at R\ = 800 and the wake behind a cylinder at R\ = 140. In 



both cases data have been recorded at about 25 integral scales downstream |1J] . In 
fig. la S'e/r^^^^ computed for the two experiments, is plotted as a function of r. In 
fig. lb we show S^/r'^^^^ as a function of r. In both figures a scaling region is observed 
only for the highest Rx. In contrast if the relative scaling is used, see fig. 2, a 
clear scaling is present for both Rx with /5(6,3) = 1.78. The vertical dashed lines 
in the fig. 2 correspond to r = brjk and they roughly indicate the extension of the 
scaling (^, that is Sr/^ < r < L. 

The ESS scaling has been checked both on numerical data and in experiments, 
in a range 30 < i?A < 2000. A direct consequence of the scaling (H) is that for all p. 
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Sp can be written in the following way: 



S,(r) = C, US 



tH-) 



(5) 



with Ul = 5*3(1/), L = Ul/e being the integral scale and Cp dimensionless constants 
selected in such a way that /(x) = 1 for x >> 1. Eq.(|^) has been carefully checked 
by computing the function 




Sp{r) \ ci) 



r \Cp US J 

If fp is independent of p this means that eq. (|^) is satisfied. 

This is seen in fig.3 where log(/6//2) and log(/4//2) are plotted as function of r/rjk. 
We clearly see that the both ratios are close to 1 within 2% for r > 5r]k. This result 
shows that eq.(^) is satisfied for 5?]^ < r < L. 

ESS has been also checked for the temperature and velocity fields in Rayleigh- 
Benard convection [l7| and in the case of a passive scalar|18|. It turns out that ESS 



is a very useful tool in order to distinguish between Kolmogorov and Bolgiano scaling 
T7|| , UTOf . In the case of the Bolgiano scaling it has been found that ({3) = 2.08 



which is clearly very different from the Kolmogorov value (^(3) = 1. In spite of this 
large difference between the values of the exponent C{^)y using ESS one discovers 
that the ratio C(p)/C(3) in the case of Bolgiano are equal to those of homogeneous 
and isotropic turbulence. The same property is observed for the ({p) obtained from 
measurements done on the solar wind for MHD. In Table I we compare the 



C{p) measured in different physical systems and the ratios 3) for MHD and 
Rayleigh-Benard convection. 

Another interesting observation concerns the behaviour of 3) with respect 
to R\. The values of P{p,3) reported in Table 1 have been measured in the range 
30 < -Ra < 5 10^. First we note that, within error bars, any change or trend of 
P{p,3) as a function of -Ra is absent. Second, we show in fig.4 the dependence of 
/3(6,3) on Rx (a similar result has been reported in ref. |2l|). This means that far 



away from boundaries the 3) are constants which do not depend on Re and on 
the way in which turbulence has been generated. 

A final point regarding ESS, concerns the generalization of the Refined Kol- 
mogorov Similarity Hypothesis. 

The RKS hypothesis states that e.^ ~ 6v^/r, as far as concern the dependence on the 
scale r, and supports eq.(|^). We can generalize the RKS hypothesis by introducing 
an effective scale L(r) = S^{r)/e, as suggested by ESS, and we obtain the following 
relation: er = 6v^ e/^s. 
Generalization of RKSH simply states that: 

Sp{r) =< MrYl >= S,{rr/\ (6) 
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In section 6 we give some theoretical support of eq.(|D. 

Eq.(^ has been first proposed in ref. [H] and carefully checked in ref. [^. A typical 
experimental result is shown in fig. 5 where < > is plotted as a function of 
5*6 (r). The energy dissipation has been computed using the 1-dimensional surrogate 
that is: 

In fig. 5 one can see a clear scaling extending over almost ten decades from the 
integral scale to rjk- The slope of the straight line is 1.005 showing that eq.(|^) is 
compatible with experimental data. 

One can argue that eq.(^) is a trivial one because for r < rj^, e,. is constant and 
Sp oc r^, thus the scaling Sn oc S^^^ is obviously satisfied. Furthermore for r in the 
inertial range eq.(§) is certainly verified because (S^/e) oc r. However in principle 
the proportionality constant of eq.(^) in the inertial and in the dissipative range 
could be different. The fact that experimentally they are found equal has several 
important consequences which will be discussed in section 5. 



3 Systems where ESS is not observed 

In the previous section we have discussed several systems where not only the ESS 
works but also the exponents /5(n, 3) are universal because they do not depend on 
the systems and on Re. We want to stress that this kind of universality, observed 
in different flows, disappears if the system is influenced by the presence of a strong 
mean shear. In this case ESS does not work, because an extended range of scaling 
is not present when Sn is drawn function of S^. 

Violation of ESS has been observed experimentally in boundary layer turbulence 



T5[| and in the shear behind a cylinder p4| . 

In a recent numerical simulation the effect of the shear on scaling laws has been 
carefully investigated using a Kolmogorov flow. 

This simulation concerns a 3D fluid occupying a volume of = sites with L = 
160, and forced such that the stationary solution has a non-zero spatial dependent 
mean velocity < v{ x) >= xsin(^), where x is the versor in the direction x, and 
L is the integral scale. 

In figures 6a) and 6b) we show the standard ESS analysis by plotting log(S'6) as a 
function of log(5'3) for two specific levels Za and Zb, where Za and Zb were chosen at 
the level of minimum and maximum shear respectively. The Rx of the simulation 
was 40 and no scaling laws were present if examined as a function of r. 
Nevertheless, it is clear from figure 6a) that ESS is observed for the case of minimum 
shear and it is not observed for the case of maximum shear (figure 6b). In both 
figures, the dashed lines are the best fit done in the range between the 20-th and 
30-th grid point and correspond to the slopes /3(6,3) = 1.78 and /?(6,3) = 1.43 for 
the minimum and maximum shear respectively. 
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However one finds that generalized Kolmogorov similarity hypothesis eq.(|^) is 
satisfied also for values of r where ESS is no longer satisfied. In order to highlight 
the previous comment we consider again the above mentioned Kolmogorov fiow. 
In fig. 7a and 7b we show the result of the scaling obtained by using eq.(^ at the 
correspondent z- levels of figures 6a and 6b for p=6. As one can see the generalized 
Kolmogorov similarity hypothesis is well satisfied in both cases although for Zb ESS 
is not observed. This is another important experimental and numerical result which 
we will consider again in section 5. 

The other relationship which we have observed to hold from large to small scales 
even in absence of ESS is the moment hierarchy recently proposed in ref. [|lT| and 
rewritten in terms of velocity structure functions in ref. p6| . 



4 Hierarchy of structure functions 



In a recent letter |TT| She and Leveque have proposed an interesting theory to explain 
the anomalous scaling exponents of velocity structure functions. The theory yields 
a prediction 

C(p)=p/9 + 2(l-(2/3)^/=^) 



which is in very good agreement with available experimental data|Tj . 

The She Leveque model is based upon the fundamental assumption on the hi- 
erarchy of the moments, < >, of the local energy dissipation. Specifically they 
consider that: 



< > " V < e"-i > ' 



where are geometrical constants and e^°°^ = lim„^oo(^^;r^) is associated in 



ref. [|TT| with filamentary structures of the flow. On the basis of simple arguments it 
is assumed that: el°°^ oc r~^/^. The value of P predicted in ref.[0 is 2/3. Notice 
that in eq.(|^) for n=l, taking into account that < e.^ >= e is constant in r, one 
immediately finds that 

(400))(1-.) ^< 4 1 

where eq.(^) has been used. 



Equation (R), which has been experimentally tested in ref. [27|, can be extended 



to the velocity structure functions|2^. Taking in eq.(|) the value n = p/3 and using 
equations (^) and after some algebra one finds the following relation for the 
velocity structure functions: 

F,^,{r) = C, {F,{r) f ■ F{r) (10) 

where 
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and 

(1-^')- 



^^3 



{1+/3) 



Cp are geometry dependent constants and = 



1/3 



Notice that eq.(|TOD is certainly valid for any /5 in the dissipative range where 
Sn oc r". Equation (^) has been experimentally tested in ref. [^. 

This can be seen in figures 8a) and 8b) where the scaling obtained for various p 
using eq. ([lO|) is reported for two different Re. As we have already observed in the 
case of eq.(P) (figure 5), the scaling extends from large to small scales even for values 
of r where ESS is no longer satisfied. 



5 A generalized form of ESS 

In sections 3) and 4) we have shown that the GKRS eq.(|^) and the hierarchy of 
moments eq.(0) are two relations which are satisfied even in fiows where ESS is not 
observed. These results suggest that the concept of ESS could be generalized in such 
a way to take into account the scaling relations equations (|^) and (p!OD properly. 
For this purpose we introduce the dimensionless structure function 



G (r) = ^^^^^ (11) 



According to Kolmogorov theory (|TT]) should be a constant both in the inertial and 
in the dissipative range, although the two constant are not necessarily the same. 
Because of the presence of anomalous scaling Gp{r) are no longer constants and by 
using @ we have: 

Gp{r) =< e^/=^ > (12) 

Thus the functions Gp{r) satisfy the hierarchies (H) and (0). Following the results 
of sections 3) and 4) equation ([T^) is valid for all scales even in cases where ESS is 
not verified. Therefore, it seems reasonable to study the self scaling properties of 
Gp{r) or, equivalently, the self-scaling properties of the energy dissipation averaged 
on an interval of size r: 

Gp(r) = G,(r)^(P'5) (13) 

where we have by definition: 

p{p, q) is given by the ratio between deviation from the K41 scaling. It will play an 
essential role in our understanding of energy cascade. Indeed, it is easy to realise 
that it is the only quantity that can stay constant along all the cascade process: from 
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the integral to the sub-viscous scales. It is reasonable to imagine that the velocity 
field becomes laminar in the sub-viscous range, Sp{r) oc r^, still preserving some 
intermittent degree parametrized by the ratio between corrections to K41 theory. In 
order to check the validity of eq.(|T^) we have plotted in fig.9 Gq^t) versus G^i^r) for 
many different experimental set-up |^,|2^, |]17|, |2^, done at different Reynolds 



numbers and for some direct numerical simulation with and without large scale 
shear. As one can see the straight line behaviour is very well supported within 
experimental errors (of the order of 3%) no deviations from the scaling regime are 
detected. Similar results are obtained, using different Gp{r) and Gq{r). 

There are two alternative ways to check ([T3|) . First of all, one can rewrite it in 
the following way: 

Spir) = (5, (r) (53 (r))P/=^-''(P'«)'?/=^ (15) 

If ([T3|) is true then p{p, q) should be equal to r(p, q). One can use (|15]) directly and 
perform a two variable fit of p(p, q) and r(p, q). Then the quantity: 



p(p, q) - r(p , q) 



= -Jz—:\ (16) 



gives a measure of the accuracy of ([T3|). We have computed Op^q for p and q in the 



range [1,8] for all the experimental and numerical results. We have found that: 

sup[apj = 0.01 (17) 



Where the above test has been done over all the experimental and numerical data 
available to us. This result tells us that the accuracy of (|1^) is extremely well 
verified. 

A second and independent check of (|13l) can be done by using (|T0|) . Indeed (|10|) 
can be checked either for fixed function of r or for fixed r as a function of 

p. In the second case we may assume that the constant Cp in (|T0|) is p-independent 
and, by plotting in a log-log scale Fp+i against Fp for fixed r and different we 
can estimate the exponent If (|TEp is true than we should observe scaling (|TD|) 



both at large scale and very small scale with the same value of Let us remark 
that the previous statement (on which our test of (|13|) is based) depends on the two 
assumptions that the log-poisson hierarchy for structure functions is true and that 
the constant Cp in (|10D are p-independent. In fig. 10a and 10b we show a log-log 
plot of -Fp+i against Fp for p = 1, ...6 and r = Sr/t and r = SOr/t for the case of two 
numerical simulations, namely Rayleigh Bernard thermal convection and channel 
flow. As one can see a clear scaling is observed with the same scaling exponent 
both for small and relatively large values of r. This confirms the quality of the 
generalized ESS scaling (^). 
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6 A theoretical interpretation 



The aim of this section is to discuss a possible theoretical interpretation of the 
experimental and numerical results previously shown. Our starting point is to revise 
the concept of scaling in fully developed turbulence. 

Let us consider three length scales ri > r2 > and our basic variables to 
describe the statistical properties of turbulence, namely the velocity difference Sv{ri). 
We shall restrict ourself to those statistical models of turbulence based on random 
multiplier. 

Thus we shall assume that there exists a statistical equivalence of the form: 

6v{ri) = aij6v{rj) (18) 

where < rj and random number with a prescribed probability distribution 

P 

By definition, we have: 

ai3 = 012023 (19) 

Equation (|T^ is true no matter which is the ratio ^ and ^. Now we ask ourselves 
the following question: what is the probability distribution Pij which is functionally 
invariant under the transformation (|l9|) ? This question can be answered by noting 
that equation (|T9D is equivalent to write: 

log ai3 = log ai2 + log 023 ( 20) 

(we assume aij > 0). Thus our question is equivalent to ask what are the proba- 
bility distribution stable under convolution. For independently distributed random 



variables a solution of this problem can be given in a complete form [^, |^. If 
the variables are correlated the situation becomes much more difficult to solve, as 
it is well known from the theory of critical phenomena. For the time being we shall 
restrict ourselves to independent random variables. 

In this case, for instance, the gaussian and the Poisson distributions are well 
known examples of probability distribution stable under convolution. These two 
examples correspond to two turbulence model proposed in literature, namely the 
log- normal model and the log- Poisson model |n|, |3l|, ||3^, |3^. A more general 
description can be found in 



We can have a different point of view on our question which is fully equivalent 
to the above discussion. A simple solution to our question is given by all probability 



distribution Pij such that: 



<•) - n ^ (21) 



for any functions gk{ri) and 7fc(p) ((■■•) represents average over Pij). Indeed we 
have: 
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We want to remark that equation (pT]) represents the most general solution to our 
problem, independent of the scale ratio Vi/rj. 

Let us give a simple example in order to link equation ( pT)) to the case of probability 
distribution stable under convolution. Following and |^ let us consider 

the case of a random log-Poisson multiplicative process, namely: 

ai, = A^(3^ (23) 

where x is a Poisson process P{x = N) = -^^^ — . 
By using (^) we obtain: 

H) = ^f,exp(a,(/5^-l)) (24) 



Equation (^J) is precisely of the form ( pTD if we write 



A, = expC., = i?2M (25) 

9i{ri) 92{rj) 

In order to recover the standard form of She-Leveque model we need to assume that 
(see also |3lD : 

giin) ~ r> (26) 

92{r,) ~ (27) 

This example highlights one important point in our discussion, i.e. the general 
requirement of scale invariant random multiplier (0) does not necessary imply a 
simple power law scaling as expressed by the equations (p^-p7|). Moreover, the 



general expression (|2lD is compatible only to infinitively divisible distribution. For 
instance, previous random multiplier model for turbulence, such as the /5-random 
model or the p-model, cannot be expressed in the general form (pTf ) independently 
of the ratio ri/rj. 

It is worthwhile to review the multifractal language at light of the previous 
discussion. In the multifractal language for turbulence, the two basic assumptions 
are: 

I) The velocity difference on scale r shows local scaling law with exponent h, i.e. 

6v{r) ~ r^; 

II) the probability distribution to observe the scaling 6v{r) ~ is given by 

In the multifractal language, therefore, there are two major ansatz: one con- 
cerns power law scaling of the velocity difference (assumption I) and the other one 
concerns a geometrical interpretation (the fractal dimension D{h)) of the probabil- 
ity distribution to observe a local scaling with exponent h. How is it possible to 
generalize the multifractal language in order to take into account equation (|2lD? 

As we shall see, the theory of infinitively divisible distribution is the tool we 
need to answer the previous question. All pubhshed model of turbulence based on 
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infinitively divisible distribution are equivalent to write D{h) in the form: 



3-D{h) = dof 



h — h 



do 



where do and ho are two free parameters while the function f{x) depends only on 
the choice of the probability distribution. For instance for log-normal distribution 
f{x) = x^. Equation (^) allows us to write: 



{Sv{r)P) = / dnihyPr 



„hp^3-D{h) _ ^hop+doH{p) 



where 



H{p) = mi{px + f{x)) 



(29) 



(30) 



We can see that equation (|29D is equivalent to a random multiplicative process given 
by: 



A [A (31) 

Equation (plf ) can be generalized to the form (^) by allowing ho and do to depend 
on r, i.e. 



ho(rj) 



Jo(rj)' 

J_ 

jioirj) 



Hip) 



where: 



hoir) = hoSh{r) do{r) = doSd{r) 
Equation (|3^ is equivalent to (^) by using: 



(32) 



(33) 



(34) 

(35) 
(36) 



(37) 



Note that the saddle point evaluation of ( pUj ) is not spoiled by the dependence of ho 
and do on r. 

We have seen that (^l]) can be reformulated in terms of multifractal language 
for infinitively divisible distribution whose function D{h) can be rewritten as in 



71 (P) =P 
72 (p) = H{p) 

The same results can be obtained by (^9]), i.e. we have 

{6v{ry) = r'^o{r)p+dQir)H(j)) 



I). We can ask the following question: what is the physical meaning of (plf ) or its 



multifractal analogous (|32|-|37|)? It is precisely the multifractal language which allows 
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us to answer this question. Indeed, the two basic assumption for the multifractal 
language can now be replaced in the following way: 

I) the velocity difference on scale r behaves as 

Hr) gi{r)g2{rr; (38) 

II) the probability distribution to observe I is g2{r)-^^^\ 
Then we have 

{SvirY) = [ dfiix)g,irrg2irY^+f^^^ = gi{rYg2{r)''^''^ (39) 



by employing a saddle point integration. The most clear physical interpretation of 
(p9|) is that the probability to observe a given fluctuation of the velocity difference 
has no more geometrical interpretation linked to the fractal dimension D{h). The 
probability distributions are controlled by a dynamical variable g2{r) which at this 
stage we still need to understand. An insight on the dynamical meaning of g2{r) 
can be obtained by the following considerations. 

Let us define e(r) the average of the energy dissipation on a scale r. We can 
define the eddy turnover time r(r) on scale r as: 

^ - «(0 (40) 

We have seen that all experimental and numerical data suggest that the following 
relation is always ( see also eq.6): 

e(r) 6v^{r) 

(e) (Sv^r)) ^ > 



where =^ means that all moments on the r.h.s. are equal to l.h.s. By using ( ^Q - ^lD 
we obtain the definition of length L(r): 

(6vHr)) 

L{r) = 6v{r)r{r) = ^ ^ " (42) 

L(r) cannot be regarded as a real length scale in the physical space. Rather, L(r) 
should be considered as a dynamical variable entering into the statistical descrip- 
tion of turbulence. This is precisely the idea behind ESS which reformulate the 
scaling properties of turbulence in terms of L{r). Indeed in order to obtain ESS 



from ( p9D it is sufficient to state that, within the range of scales where ESS is ob- 
served, (yfi(r)^/'^" ~ g2{rY^^^^ ~ L{r). The physical meaning of ESS is strictly linked 
to (^) and in particular to (^Tj) which is a generalization of Kolmogorov Refined 
Similarity Hypothesis. 

Let us summarize all our previous findings: 

A) we have introduced the idea of scale invariant random multiplier satisfying equa- 
tion (^; 
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B) we have shown that infinitively divisible distributions are all compatible with(pTD; 

C) we have shown that the multifractal language specialized for the case of infini- 
tively divisible distribution gives equation (pT]) (with n = 2 and 71 (p) linear in p) 
and it is equivalent to scale invariant random multiplier; 

D) finally we have argued that the correct scaling parameter to describe the sta- 
tistical properties of small scale turbulent flows is not directly linked to a simple 
geometrical interpretation, rather it should be considered a dynamical variable. 
Our finding A-D enables us to have a unified theoretical interpretation of the ex- 
perimental and numerical results presented at the beginning of this paper. Indeed 
equation (|37| ) or (^) tells us that the anomalous part of the structure functions: 



'Gp\n = /r...^/^AVW.^ (43) 



satisfies the scaling properties: 



Gp(r) = G,{rY^- (44) 



where pp^q = ^z^- According to our analysis of the experimental and numerical 
results, the scaling (44) is observed down to the smallest resolved scale. 
We have shown that, in the theoretical framework so far exposed, we recover the 
ESS when ^i(r)i/'^« ~ gl^'^° ~ L(r). If ^i(r)i/'*o ^ ^2(r)^/'^» we lose ESS, but its 
generalized version (H3) is still valid. 



6.1 Synthetic Turbulence 

We can also use (0)and (|39| ) to simulate a synthetic signal according to a random 
multiplicative process satisfying (|32D . This can be done by using the algorithm 
recently introduced in |^ 



Let us consider a wavelet decomposition of the function 0(x): 

00 

0(x) = ^ aj-fc^j-fc(x) (45) 

j,k=0 

where ipj^ki^) = 2^^'^ip{2^x — k) and ip{x) is any wavelet with zero mean. The above 
decomposition defines the signal as a diadic superposition of basic fluctuations with 
different characteristic widths (controlled by the index j) and centered in different 
spatial points (controlled by the index k). For functions defined on = 2" points 
in in the interval [0, 1] the sums in (^) are restricted from zero to n — 1 for the 
index j and from zero to 2-' — 1 for k 



In [0 it has been shown that the statistical behavior of signal increments: 

< |50(r)|P >=< |0(x + r) - (t){x)\P >~ r^^^^ 
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is controlled by the coefficients aj^k- By defining the a coefficients in terms of a 
multiplicative random process on the diadic tree it is possible to give an explicit 
expression for the scaling exponents C(p)- Fo^^ example, it is possible to recover the 
standard anomalous scaling by defining the a's tree in term of the realizations of a 
random variable t] with a probability distribution P{r]): 



"0,0 



«i,o = ^71,0 "0,0; = '7i,i«o,o; 

"2,0 = ^?2,0ai,0; "2,1 = ^?2,iai,0; "2,2 = ^?2,2ai,i; "2,3 = ^?2,3ai,i, (46) 

and so on. Let us note that in the previous multiplicative process different scales 
are characterized by different values of the index j, i.e. rj = 2~K If the rjj^k are i.i.d. 
random variable it is straightforward to realise that Uj^k are random variables with 
moments given by: 

< >=r-^°^^(^) (47) 

where the "mother eddy' q;o,o has been chosen equal to one. In p7|) with we 
intend averaging over the P{ri) distribution. In p3 it has been shown that also the 



signal 0(a;) has the same anomalous scaling of (0). 

In order to generalize this construction for function showing ESS or generalized- 
ESS scaling of the form (|37|) and (^4]) is now sufficient to take a probability dis- 
tribution, Pi (77), for the random multiplier with the appropriate scale dependency 
(pip. This will be implemented by allowing a dependency of P{rjjk) on the scale 
Tj = 1~\ i.e. the 77's random variables will be still independently distributed but 
not identically distributed with respect to variation of the scaling index j. 
According to the previous discussion, ESS corresponds to have only one seed- 
function defining the multiplicative process, i.e. g\{r) the range 
of scales where ESS is valid (r < br]k)- On the other hand, at scales smaller then 
5 ~ 6 Kolmogorov scale, ESS is not more valid because g^'''^ begins to deviate sub- 
stantially from 5^1^'^°: only G-ESS should be observed and we need a multiplicative 
process defined in terms of two different seed-functions. 
Following this recipe we define the signal such that : 

< 5v{rY >= U^F{ry/^G{r)^^P'>-P/''^ (48) 

where 

F{r) =< 6v{r)^ > /U^ (49) 

The function G{r) is defined in such a way that for F{r) much greater than rjk/L, 
G{r) ~ F{r) while for very small scales r we have G{r) ~ i]k/L. In the following we 
choose the simplest ansatz: 

G{r) =B + AF{r) (50) 
with B = T]k/L and A is a dimensionless constant. 
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Let us now spend some words in order to clarify the previous definitions. Relation 
(^) is defined such that experimental results are reproduced with good accuracy and 
G-ESS scaling (|4^) is satisfied by definition. By assuming (|50|) the only unknown 
function is F{r) =< Sv{r)^ > /Uq. On the other hand the function < 6v{r)^ > is 
always very well fitted by the Batchelor parameterization: 

r/3 3 

< Sv{r)^ >= r^^. (51) 

From expression (^8[) is immediate to extract the expression for the two seed- 
functions gi{r),g2{r) used in the previous sections, namely: 

^2(r) = G{rY° (52) 

Let us note that gi{r) goes smoothly from the intermittent value, gi{r) ~ {ho = 
1/9 for the case of She Leveque model), assumed in the inertial range to the laminar 
value, (71 (r) =~ r, characteristic of scales much smaller than Kolmogorov scale. 

For the practical point of view we have constructed our signal by using a random 
process for the multiplier rjj,k{i^j) with a scale-dependent Log-Poisson distribution. 
The scale dependency of parameters entering the distribution has been fixed in 
terms of relations (^2|) and and such that the ({p) exponents correspond to the 
She-Leveque [|ll| expression, namely: 



vAr,) = (53) 

where Poisson variable with mean Cjj+i = log{g2{rj+i) /g2{rj)), Ajjj^i = 

gi{rj)/gi{rj+i) and /? = 2/3. This choice leads to the standard Log-Poisson scaling 
in the inertial range: 

(1 - (3P/^) 

C(p) = V+(1-3M 

and to the following expression for the ratios of deviations to the Kolmogorov law: 

H{p)-p/3 1-/3-/^ 
P^'^- H{q)-q/3' ""^^^ (1-/3) 

Signal constructed according to this scenario will be referred to as signal- A in 
the following. 

In fig. 11 we show the structure function of order 6 for such a signal plotted versus 
the separation scale r at moderate Reynolds number. Clearly, for this choice of 
Reynolds number there is not any inertial range of scale where scaling exponents 
could be safely measured. On the other hands, our signal shows G-ESS scaling, how 
it is possible to see in fig. 12. 
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7 Multiscaling 



We now turn our attention to a different question wliicli is connected to tlie the- 
oretical results so far discussed, namely the role played by viscous effects. It is 
generally argued that the anomalous scaling can be observed for scale larger than a 
given viscous cutoff. The physical interpretation of this statement is that non linear, 
intermittent, transfer of energy is acting only for scale larger than the viscous scale. 
Below such a scale the structure function are supposed to show a simple (regular) 
scaling {Sv^{r)) ~ r^. 

Usually the viscous cutoff is introduced as the scale at which the local Reynolds 
number is of order one, namely: 

^ ~ 1 (54) 

This condition can be obtained by the requirement that the local energy transfer 
e(r) ~ (5f^(r)/r) becomes equal to the energy dissipation z/(5f ^(r)/r^): 

5f^(r) 6v^(r) 
V ~ — - — (55) 



which gives equation (|54|) . There is a well defined prediction, based on (pif ), for- 
mulated by Frisch and Vergassola using the multifractal language. Indeed for any 
exponent h one can introduce the /i-dependent viscous cutoff given by: 

- y (56) 

where (5f (r) ~ . 

It follows that Tdiji) is a fluctuating quantity. There are two consequences of this 
theory. The first one predicts that for the structure functions {6v^{r)) there exists 
a cutoff scale Vp dependent on p and moreover Vp < Vg for p > q. 
The second prediction concerns the moment of the velocity gradients F which are: 

(FP) ~ {rdihY'^-^Kdihf-'^^^^) ~ Re-'^P^ (57) 

with z{p) = snpj, [{{h - l)j9 + 3 - D{h)) /(I + h)]. 

Between the two predictions the first one is qualitatively more peculiar of (pU]). 
In particular, the first prediction states that between the end of the inertial range 
(i.e. the region where anomalous scaling of {5vP{r)) with respect to r is detected) 
and the dissipation cutoff r^, the local slope is controlled in rather complicate way by 
D{h). The second prediction is somehow weaker because present experimental data 
do not distinguish among several models, so far proposed, for the Re-dependence of 

In order to compare the multiscaling in the dissipation range with our experi- 
mental and numerical data, we have produced a synthetic turbulence signal (signal 
B hereafter) similar to the one already discussed but with and independent 
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on r. The effect of dissipation is introduced by using (|5J). In fig. 13 we compare 
the local scaling exponents d{[og{5v^)) / d{\ogr) for p = 6 between the two synthetic 
signals. 

In fig. 14 we plot the relation (^) for p = 6 for signal A and B. Finally in fig. 15 
we plot the ratio G^jG^^''^ for signal B. By comparing figures 13, 14 and 15 with 
the analogous experimental and numerical results discussed in the previous sections 
(see fig 5,7,9), we can state that the quantitative and qualitative prediction based 
on ( p6D is not verified by experimental data. On the other hand, signal A, based 
on an explicit r-dependence of and /iq, seems to be more closely related to what 
observed experimentally. Let us remark that signal A has no cutoff effect imposed 
by condition ([5^). 

The above discussion rules out the effect of multiscaling on the viscous cutoff 



Previous claims on the validity of multiscaling effects should be considered either 
wrong or affected by experimental errors. On the other hand our model, used to 
implement the synthetic signal A, should be considered a very accurate model even 
for scale close to the regular region where bvir) ~ r. 

There is, however, a theoretical question concerning multiscaling which we are 
still not able to answer completely and that we shall try to formulate in the following. 
There are two possible scenario in which a viscous cutoff may be considered. 

In the first scenario (let us call it scenario I) we can imagine to consider equation 
(0) as a fundamental relation independent of any other theoretical considerations. 
The idea is that when the local Reynolds number is sufficiently small, then non 
linear effect must be neglected. In order to compute the viscous scale, one should 
make use of the relation: 

bv{r) ~ \J^F{rYI''G{rf-^l'' (58) 

obtained by the two definitions (|38| ) and (^). In the equation (|^) h is now the 
standard multifractal scale-independent exponent. Generalized scaling (^) should 
be considered realized with probability 

P,(r) ~ G(r)3-^W (59) 



inserting (|58D into (|5J) we obtain: 

r^\J^F{T^fl^G{r^f-^l^ ~ v (60) 
where Vdih) is the fluctuating cutoff. 

We now look for a solution of equation (^) in a region where 



^3 



F(r) ~ (61) 
G{r) ~ I (62) 
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After some algebra we obtain 



3/1-1 



r,(/i) ~ Re—7]k (63) 

Thus rd{h) is a fluctuating quantity as in (^). These fluctuations, however, happen 
in the region where Sv{r) ~ r and therefore no effect on the scaling of structure 
function is produced. From (58-^) we can compute the scaling of (F^*) as function 



of Re. The scaling is independent on rd{h) and it is: 

(TP) ~ i?et(P-^'') (64) 



consistent with ( ^8l) in the limit r 0. 

Note that (F^) ~ i?et(2-C2)_ xhis implies that for (2 7^ 2/3, (F^) does not scale as 
Re. If we want to recover the experimental fact that (e) is constant with Re, we 
should allow for a Re-dependent constant in (pHj). 

At any rate, because C2 — 2/3 is a small quantity, these effects are quite small in the 
full range of available Re-number. 

We can summarize the scenario I as follows: the scaling (|4^) and ( ^81) are verified to 
all scales; the condition ( |5^ ) introduces a viscous cutoff which fluctuates in the region 
where 6v{r) ~ r; intermittency in the gradient of the velocity field are prescribed 
by (0). 

The above conclusions imply that scaling (^) must be violated near the viscous 
cutoff, as one can immediately check by an explicit computation. One can take an 
opposite point of view and assume that (0) is a fundamental relationship which 
must not be violated. This corresponds to the second scenario. 

In the second scenario (E^ is disregarded and one generalize (E^j as: 



^r^ ~ (65) 

where F is the velocity gradient and is the viscous cutoff. In order to compute 
Td, one observes that F ~ where r^Vd) is the eddy turnover time at the viscous 
cutoff. We obtain: 

^ ^ Sv^{rd)T{rdy ^ 5v{rd){5v^{r-d)) ^^^^ 
{Sv^{rd)) e 

where following (^), we used T{rd)Sv{rd) = {6v'^{rd)) /e. Once again, by using (^8]), 
(|6l|) and (|62|) we can obtain an explicit formula for rd{h): 

rd{h) ~ Re-^^^"^Su (67) 

Thus also in scenario II, we have strong fluctuations of the viscous cutoff. 
The computation of the gradients is quite straightforward from (p5|). We have: 
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where we have used (|5^). Finally by using (^) and (^2]) we get: 



-p/2) 



(69) 



Note that in scenario II (F^) ~ Re because ^3 = 1- 

Because of (|69[) , the II scenario violates (|58|) and for scales smaller than rrf(/i) 
while ( Pip is always satisfied. 



It is quite difficult to understand which one of the two scenario is actually verified 
by experimental and numerical data. In most cases the scale resolution does not 



reach the region where 6v{r) ~ r. At any rate, either (^) or p^) should be violated 
at very small scales as the result of viscous effects. This violation is rather small 
and may not be easily detectable at low or moderate Reynolds numbers. 
The common point about the two scenario is that the viscous cutoff (if any) acts 
at scale where already the velocity structure functions behave in a regular way, i.e. 
5v{r) ~ r. 



8 Conclusions 

In this paper we have proposed several new results concerning the scaling behaviour 
of small scale statistical properties of turbulence. It is worth to summarize our main 
findings trying to outline questions which are still to be answered. 

1) We have reviewed the main results on ESS and in particular we have shown 
that in homogeneous and isotropic flows in turbulence, Rayleigh Benard convection 
and solar wind magnetohydrodynamics, the ratio C(p)/C(3) seems to have an uni- 
versal behaviour. This is a rather striking and unexpected result which implies that 
anomalous violation of dimensional scaling may be explained in an universal way. 
We do not know any simple phenomenological explanation for our finding. 

2) We have shown that ESS is not observed when relatively strong shear flows are 
present. A phenomenological analysis, based on the Kolmogorov equation, shows 
the relevance of a length scale based on the mean energy dissipation and the shear 
strength. This analysis should be refined in order to acquire more quantitative 
predictions. At any rate, our observation suggests that previous finding of violations 
of ESS should be due to the presence of shear flows. 

3) We have shown that the refined Kolmogorov similarity can be generalized by 
including ESS. This generalization is verified extremely well in both experiments 
and numerical simulations. More important, we have shown that the generalized 
refined Kolmogorov similarity is true also in cases where ESS is not observed. 

4) Similar to the previous point, we have shown that the hierarchy relation 
based on log-poisson distribution for the structure functions is very well supported 
by experimental data, also for very small scales where ESS is not observed. 
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5) Based on our results in l)-4) we have proposed a generalization of ESS. This 
generalization is supported both by experimental and numerical data and it seems 
not affected by viscous cutoff. 

6) We have developed a theory which unifies the previous point. The theory 
is based on the assumption that the probability distribution is infinitively divisible 
and predicts the existence of the generalized ESS. The theory can also be used to 
generate artificial signals which displays all the scaling features observed in real 
data. 

7) We have shown that the original proposal on the multiscahng for the viscous 
cutoff is incompatible with the turbulence data. The theory formulated in this paper 
removes this incompatibility and suggests that multiscahng is acting at much smaller 
scales than previously proposed. The new point on the theory is a change of view in 
the probability distribution of the original multifractal model which is not directly 
linked to a geometrical interpretation in terms of fractal dimensions. 

8) Finally we have shown that violations of either the generahzed refined Kol- 
mogorov similarity or the generalized ESS should occur at very small scale. Our 
present data analysis does not allow us to distinguish among the two possibilities. 

9 Acknowledgements 

We thank B. Castaign, S. Fauve, U. Frisch, Z.S. She, K.R. Sreenivasan, G. 
Stolovitzky, , for many useful discussions. During the writing of this paper, we 
became aware of a manuscript of C. Meneveau dealing with a possible explanation 
of ESS. His conclusions are physically close to part of the discussion in section 6. 

This work has been partially supported by EEC contracts ERBCHBICT941034, 
CT93-EVSV-0259 and CHRX-CT94-0546. 

References 

[1] A.S. Monin, A.M. Yaglom, " Statistical Fluid Mechanics " (MIT Press, Cam- 
bridge, MA, 1975). 

[2] A. N. Kolmogorov, C. R. Acad. SCI USSR 30, 299 (1941) 

[3] F. Anselmet, Y. Gagne, E. J. Hopfinger, R.A. Antonia, J. Fluid Mech. 140, 63 

(1984) 

[4] A. Vincent, M. Meneguzzi, J. Fluid Mech. 225, 1 (1991). 

[5] A. N. Kolmogorov, J. Fluid Mech. 13,83 (1962). 

[6] U. Frish, R Sulem, M. Nelkin, J. Fluid Mech. 87, 719 (1978) 



20 



[7] G. Parisi, U. Frish, " On the singularity structure of fully developed turbulence 
" , in: "Turbulence and predictability in geophysical fluid dynamics and climate 
dynamics", ed. M. Ghil, R. Benzi, G. Parisi(North Holland, Amsterdam, 1985), 
p. 84 

[8] R. Benzi, G. Paladin, G. Parisi, A. Vulpiani, J. Phys. A, 17, 3521 (1984). 

[9] C. Meneveau and K. R. Sreenivasan, Phys. Rev. Lett., 59, 1424, 1987. 
[10] Y. Gagne and B. Castaing, C.R. Acad. Sci. Paris, 312, 441, 1991. 
[11] Z. S. She, E. Leveque, Phys. Rev. Lett. 72, 336 (1994). 
[12] U. Prish, M. Vergassola, Europhys. Lett. 14, 439 (1991) 

[13] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, S. Succi, Phys. Rev. E 48 
(1993) R29. 

[14] R. Benzi, S. Cihberto, C. Baudet, G. Ruiz Chavarria, Physica D 80, 385 (1995). 

[15] G. Stolovitzky, K. R. Sreenivasan, Phys. Rev. E, 48, 32 (1993). 

[16] M. Briscolini, P. Santangelo, S. Succi, R. Benzi, Phys. Rev. E 50, 1745 (1994). 

[17] R. Benzi, R. Tripiccione, F. Massaioli, S. Succi, S. Ciliberto, Europhys. Lett. 
25, 331 (1994). 

[18] G. Ruiz Chavarria, C. Baudet, S. Cihberto, " ESS in passive scalars " to be 
published in Europhys. Lett. (1995). 

[19] S. Cioni, S. Ciliberto, J. Sommeria, " Scaling laws in turbulent thermal convec- 
tion in low Prandtl number fluids " to be published in Europhys. Lett. (1995). 

[20] R. Grauer, Physics Lett. A 195, 335 (1994). 

[21] F. Belin, P.Tabeling, H. Willaime, " Exponents of structure functions in a low 
temperature helium experiment ", submitted to Physica D. 

[22] G. Ruiz Chavarria, J. Physique, 4, 1083 (1994). 

[23] G. Ruiz " Scaling law in boundary layer turbulence " in preparation (1995). 

[24] R. Benzi, S. Ciliberto, C. Baudet, G. Ruiz Chavarria, R. Tripiccione, Europhys. 
Lett. 24, 275 (1993) 

[25] R. Benzi, M. V. Struglia, R. Tripiccione " ESS in numerical simulations of 3D 
anisotropic turbulence " preprint 1995, submitted. 

[26] G. Ruiz Chavarria, C. Baudet, R. Benzi, S. Cihberto, J. Physique 5, 485 (1995). 



21 



[27] G. Ruiz Chavarria, C. Baudet, S. Ciliberto, Phys. Rev. Lett. 74, 1986 (1995). 
[28] G. Amati, R. Benzi, S. Sued, in preparation (1995). 

[29] Feller " An introduction to Probability theory and its applications ", vol.2 
Wiley, New York (1966). 

[30] Novikov E.A. Phys. Rev. E 50, R3303 (1994) 

[31] B. DubruUe Phys. Rev. Lett. 73, 959 (1994). 

[32] Z.-S. She and E.G. Waymire, Phys. Rev. Lett. 74 262 (1995) 

[33] B. Gastaing, " The temperature of turbulent flows " to be published in J. 
Physique. 

[34] R. Benzi, L. Biferale, A. Grisanti, G. Paladin, M. Vergassola and A. Vulpiani, 
Physica D 65 352 (1993). 

[35] Y.Meyer Actualite Mathematiques, Hermann Editeurs 1990. 



22 



Figure captions 

Figure 1: Structure functions Sp{r) / r'^^P'> as a function of r. Data taken from 
an experiment on a jet at -Ra = 800 (diamonds). Data taken from tlie wake beliind 
a cylinder at -Ra = 140 (triangles), (a) p=6 and ({6) = 1.78. (b) p=3 and ({3) = 1. 
Logarithms are in base 10. 

Figure 2: Structure functions Sq as a function of S3 at -Ra = 800 (a) and 
-Ra = 140 (b) computed from the same data-set of fig.l. Vertical dashed lines 
indicate the value of 5*3 at 5//^. 

Figure 3: Logarithm of ratio of the universal functions /n//2 for two cases 
n = Q (diamonds) and n = 4 (circles) for the wake behind the cylinder of fig.l as a 
function of r/rjk. 

Figure 4: Dependence of the exponent /?(6,3) as a function of Re. (-Ra — 
1.4-Re^'^^). The last point is from ref. See also ref. [pT| . 

Figure 5: Log- log plot of < > 6*3 (r)^ against 5*6 (r) at -Ra = 500. The 
straight line refers to the slope 1.005. Data are from an experiment of turbulence 
behind a cylinder and the measurement point was at about 25 diameter down stream. 

Figure 6: a) Log-log plot of ESS scaling for the longitudinal structure function 
5*6 (r) versus 5*3 (r). Data are taken from a numerical simulation of a shear flow at 
-Ra = 40. The dashed line is the best fit with slope 1.79. Every point in the plot 
corresponds to a grid point and the lattice spacing is ~ Irjk wide. The computation 
of the structure functions is performed in points of the flow where the shear has a 
minimum. 

b) the same of (a) but for points where the shear is maximum. The dashed line 
is the best fit with slope 1.43. At variance with previous case ESS is not observed. 

Figure 7: Check of eq.(^ for p=6 using the same numerical simulation of the 
shear flow discussed in fig. 6 (log-log plot). Energy dissipation has been computed 
by using the 1-dimensional surrogate in order to compare this result with laboratory 
experiments (see fig. 5). (a) points of minimum shear, (b) points of maximum shear. 
The points refer to the scales at 2, 4, 5, 8, 10, 16, 20, 32, 40 grid points and the dashed 
line is the best fit done over these points, corresponding to the slope 0.99 for both 
minimum and maximum shear data. Although in this case ESS is not observed (see 
fig. 6b), the generalized refined Kolmogorov hypothesis eq.(^) works within 3%. 

Figure 8: The function Fp^i{r) defined in eq.(|lO]) is plotted, for several values 
of p, as a function of [{Fp{r)f ■ F(r)] with P = 2/3, p' = (3^ and 5=1/3. 
-Ra = 140 in (a) and -Ra = 800 in (b). The 5 curves in (a) and (b) correspond to 
p = 1,2,3,4,5 starting from the bottom lines. They have been vertically shifted 
of —0.4, —0.2, 0, 0.2, 0.4 in order to separate them. The solid lines have slope 1. 
Logarithms are base 10. 

Figure 9: Log-Log plot of GqIt) versus 6*5 (r) for different laboratory and 
numerical experiments. (-I-) Data are taken in a wake behind a cylinder where 
standard ESS was not observed [^. (o) Data taken from the region with log-profile 
of a boundary layer (courtesy of G. Ruiz Chavarria) where standard ESS was not 



23 



observed, (squares) Data taken from a numerical simulation of thermal convection 
Tn where standard ESS was observed, f A ) Data taken from a direct numerical 



simulation of a channel flow where standard ESS was not observed 28^ 



Figure 10: Log-log plot of Fp+i against Fp for p = 1, ...6 and r = 37]^ (squares) 
and r = SOrjk (circles) for the case of two numerical simulations, namely Rayleigh 
Bernard thermal convection (a) and channel flow (b). As one can see a clear scaling 
is observed with the same scaling exponent (3' both for small and relatively large 
values of r. This confirms the quality of the generalized ESS scaling ([T3|) . 

Figure 11: Log-log plot of the 6th order structure function for the signal A 
with 19 fragmentation at small Reynolds number. Notice the absence of any scaling 
range. 

Figure 12: G-ESS (log(G6(r)) vs log(G'4(r))) for the signal A at the same 
Reynolds number of figure 10. The slope is p^^Q = 0.241 in perfect agreement with 
the theoretical prediction obtained from {1^) where for ({p) we have used the She- 
Leveque expression. 

Figure 13: 6th order local scaling exponents for the signal A (circles) and signal 
B (squares). Notice that the qualitative behaviour of the two signals is almost the 
same: both of them go from an intermittent scaling (C(6) ~ 1.8) at large scale to a 
laminar scaling (C(6) = 6) at small scales. 

Figure 14: Generalized-Kolmogorov refined hypothesis ( ^T] ) for the 6th order 
structure function in both signal A (squares) and signal B (circles). Notice the 
sudden jump at the Kolmogorov scale present when multiscaling is valid (signal B). 

Figure 15: Compensated slope (G'4(r)/G'g'''''(r)) for signal B. Notice deviations 
of the order of 10% while for signal A the same quantity is constant by definition. 
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Table I 





1^ 




(,(Bolg.) 


(,(MHD) 


J(p.3)(Bolg) 


/iOA3)(MHD) 




1 


0.37 


0.77 


0.28 


0.37 


0.36 




2 


0.70 


1.46 


0.55 


0.70 


0.70 




3 


1 


2.08 


0.78 


1 


1 




4 


1.28 


2.66 


1 


1.28 


1.28 




5 


1.54 


3.20 


1.20 


1.54 


1.54 




6 


1.78 


3.70 


1.39 


1.78 


1.78 




7 


2.00 


4.16 


1.58 


2.02 


2.02 




8 


2.23 


4.63 


1.75 


2.24 


2.24 



Table 1: We show some measured values of ({p) and 3) for 1 < p < 8. In 

the second cohimn we report the ((p) measured in 3D homogeneous and isotropic 
turbulence(30 < R\ < 2000), in the third column the ({p) measured in Rayleigh 
Benard convection when the Bolgiano scaling is the relevant one{Rx ~ 30), in the 
fourth column the ({p) obtained from the measurements of the solar wind {Rx ~ 
5 10^). We note that the C(3) of the last two cases are clearly very different from 1 
which is the value of ({3) in the second column. The ratios 3) computed from 
the values of the third and fourth column are shown in the fifth and sixth columns 
respectively. The /3(p, 3) are equal within error bars to those of the first column. 
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